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Abstract 

Two new one-dimensional fermionic models depending on two independent 
parameters are formulated and solved exactly by the Bethe-ansatz method. 
These models connect continuously the integrable Hubbard and supersym- 
metric t-J models. 
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The Hubbard model together with the t-J model are the most studied models describing 
strongly correlated electrons. In one dimension they are paradigm of exact integrability in 
the physics of strongly correlated systems. In these models we have beyond a hopping term 
t (kinetic energy) a on-site coulomb interaction U, in the case of the Hubbard model [1], or 
a spin-spin interaction J, in the case of the t-J model [2-4]. 

An interesting question in the arena of exact integrable models, we want to solve in 
this Letter, concerns the existence of a general exact solvable model containing these two 
well know models as particular cases. After the exact solution of these models [1-4] , several 
extensions which keep exact integrability were proposed, either by introducing correlated 
hopping terms [5-11], or by including an anisotropy (q-deformation) [12-15] (see ref. [16] for 
a review). However none of these extensions contains simultaneously the Hubbard and t-J 
models as particular cases. In this Letter we present two new integrable two-parameter 
models having this nice property. These models contain as particular cases the Hubbard 
model [1] and the Essler-Korepin-Schoutens model [6] as well as its q-deformed versions 
[12,14,15]. We remind the reader that the latter model [6] contains the supersymmetric t-J 
model in a particular sector. 

Our starting point is the introduction of a general one-dimensional Hamiltonian contain- 
ing all the possible nearest-neighbour interactions appearing in different exactly integrable 
models with four degrees of freedom per site. This Hamiltonian thus contains correlated 
hopping terms in the most general form, spin-spin interactions as in the anisotropic version 
of the t-J model, as well as pair hopping terms and three- and four-body static interactions 
between electrons. The Hamiltonian is given by 

L 

H j,k = Yl ( c j> c M + h.c.)[l + t al n jf ) + t a2 n kfS + t' a n jl3 n kp ] 

a(#/9) 

+ H { Jc U c k,p c 3,P c k^ + Vapn jta n ktP ) + Un jtl n jt2 

+ tpic^c+zCk^Ck,! + /i.e.) + V^n jt2 nk,in k ,2 + vf' 'n,- in M n fc)2 
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+ Vg ( 'n jtl n jt2 n k 2 + >n jt in jt2 n Kl + V A n jtl n jt2 n k)1 n k ^ (1) 

where Cj >a and n Jja = c^ a Cj a (a = 1,2) are the standard fermionic and density operators. 
The physical relevance of such Hamiltonian is discussed, e.g. in [17,18]. 

In (1) we have included a correlated- hopping interaction in its most general form, which 
depends on t a i,t a2 and t' a (a — 1, 2). In the theory of exactly integrable systems, models with 
such kinetic terms were first studied in [5,19] and the possible physical relevance of them 
is given in [20]. In the limit t a p = —t' a = — 1 this term gives a constrained hopping term 
and the condition for integrability gives the anisotropic t-J model at J = e~^Vi 2 = e'V^ = 
±1, t p — U = V3 = V4 — Q . The Hubbard model is obtained by destroing the correlation 
in the hopping term (t aj3 — t' a — 0) and by setting t p = J = V\ 2 = V21 = V3 = V4 = 0. For 
the case where J = the conditions for integrability has been investigated in [9,10], and a 
two parameter generalization of the correlated- hopping model has been contructed in [9,11]. 
Recently some one-parameter models with J ^ has been constructed [14,15,21] on the base 
of solutions of Yang-Baxter equations of vertex models [14,22,23]. In this letter we present 
the results of our investigation on thr integrability conditions in the case J 7^ 0, V aa = 
and t a/ 3 7^ 1. 

We require the wavefunctions of the Hamiltonian (1), with n electrons, to be given by 
the Bethe Ansatz 

\n >= y(r Ql ,a Ql ; . . . ; r Qn , a n )\r Ql , ...,r Qn > 
Q 

n r 

*(n, ari; ...; r n , a n ) = ^ ^p^'.pT" II x p/ > x i = expfa'fy) (2) 

p j=i 

where Q is the permutation of the n particles such that 1 < tq 1 < tq 2 < ... < tQ n < L, and 
a = 1, 2 denotes the kind of particles (up or down spin). The sum is over all permutations 



For these parameters the number of double occupied sites are conserved and the model corre- 
sponds to that introduced in [6]. The t-J model is obtained in the sector where there is no double 
occupied sites. 
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P = [Pi...P n ] of numbers 1,2, ...,n. The coefficients Ap^'p" Qn from regions other than 
Rq = [fQi < ••• < rg„] are connected with each other by the elements of the two-particle 
S-matrix 

E (3) 

q',/3'=1,2 

As a necessary condition for integrability of the model under consideration, the two-particle 
scattering matrix has to satisfy the Yang-Baxter relations [24,25]. Although we have not 
solved this problem in the general case we were able to establish exact integrability of (1) 
in the two new cases, which we denote by models A and B: 
A) 

t\ = et 2 = t 3 = et A = sinfi; t 5 — e 



J = - £tn = --U = V 12 e 2r > = V 21 e- 2r > = costf; 



2 

v aa = y 3 (1) = = v s (3) = V w = Va = o, ( i) 



B) 



ti = et 2 = et 3 e 2r> = t 4 e 2v = sin^; t 5 — e 
J = -et p = V 12 e 2v = V 21 e- 2v = costf; 
U = 2t p + ^^-ee^) 2 
V aa = V 3 {2) = V 3 {4) = V 4 = 0, V 3 {1) = -V 3 {3) = V 12 - V 21 , (5) 

where in (4) and (5) we denote 

tn = £4 — 1; tu = t 3 — 1; t 2 \ — t\ — 1; t 22 — t 2 — 1; 
t\ — t$ — t 3 — t 4 + 1; t' 2 — t§ — t\ — t 2 + 1, 

e = ±1 and t? and 77 are free complex parameters. 

For — > both cases reduce to the anisotropic t-J model studied in [12], which is the 
generalization of the supersymmetric t-J model [2-4]. More exactly, in this limit we obtain 



the q-deformed extended Hubbard model [6]. Moreover from (1,4,5) we see that the model 
B with , T) = i$,e = +1, and the model A with r\ = 0,e = — 1, reduce to the nontrivial 
q-deformations of the extended Hubbard model considered in [14] and [15], respectively. 
These models have been constructed on the base of solution of the Yang-Baxter equation 
for the R-matrix which was found by [14,23]. In the opposite limit $ — > ix/2 both models 
with e = l give us the Hubbard model, provide in model A r\ = [ln({7') — ln(cosi?)]/2, and 
in model B r] = \^U\d - tt/2\. 

The non-vanishing elements of the two-particle S-matrix of both models satisfy 

qaa -i . qctf) q(3a 

D aa — L i D a/3 ~ D /3a > 

s %i( x u x*)S$l(x2, xi) = -S$(x 2 , x 1 )S^(x 1 ,x 2 ), (6) 

and for the different models are given by 
A) 

S™p(x 1 ,x 2 ) = (xi - x 2 )b 12 (x 1 ,x 2 )/a 1 (x 1 ,x 2 ); 

SpP(x!,x 2 ) = [c (x 1 ,x 2 ) + h (x 1 ,x 2 )x 1 + b 2 (x 1 ,x 2 )x 2 - gx 1 x 2 ]/a 1 (x 1 ,x 2 ); (7) 

B) 

S%%(x 1 ,x 2 ) = {xi - x 2 )bi 2 {xi,x 2 )/a 2 {xi,x 2 )] 

SpP(x!,x 2 ) = [c (x 1 ,x 2 ) + {x x e~ 2l] + x 2 e 2r ')b 12 (x 1 ,x 2 )]/a 2 (x 1 ,x 2 ); (8) 
where a < f3 and 

ai{xi,x 2 ) = c (xi,x 2 ) + [bi(x!,x 2 ) + b 2 (x 1 ,x 2 )]x 2 - gx\\ 
a 2 (x 1 ,x 2 ) = c (x 1 , x 2 ) + (e 2v + e~ 2v )bi 2 (xi,x 2 )x 2 ; 
h( Xl ,x 2 ) = {tl + eJ 2 e~ 2r ')D 12 + Je- 2 "(xi+x 2 ); 
b 2 ( Xl ,x 2 ) = {t\ + eJ 2 e 2v )D 12 + Je 2 \ Xl + x 2 ); 
h 2 (x 1 ,x 2 ) = eDyi + J(x 1 + x 2 )\ 
Co(xi,x 2 ) = (U - 2t p )x 1 x 2 + [t p D 12 - Xi- x 2 ]D 12 ; 
D 12 = 1 + xxx 2 ; g = cos d sin 2 ^(e - ee~ r ') 2 . 
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a=l 

m 



nsn(^)= n Z ?d (i) ? ^' = 1 '-'-' (9) 



To complete the proof of the Bethe anzatz (2) we must check the eigenvalue equations in 
the sector where the total number of particles is n = 3, 4. This gives a complicated system 
of equations. A manipulation of this problem on a computer gives us the values of the 
coupling constants and V 4 in eqs. (4) and (5). The periodic boundary conditions on the 
lattice with L sites lead us to the Bethe-ansatz equations. In order to obtain these equations 
we must diagonalize the transfer matrix of a related inhomogeneous six-vertex model with 
Boltzmann weights (6-8). This latter problem can be solved by standard algebraic methods 
[26]. The Bethe-ansatz equations are written in terms of the variables Xj (xj = exp(ikj)) 
and additional spin variables 

For both models we have 

m 

(x J -) L = (-l) B - i n^i,4 1) ); J = l,-,n; 

J '- L -<■'<• J- A Ql2/„(1) 

where m < L is the number of particles with up spins. The eigenenergies of the system are 
given by 

n 

E=-'£(x j + xJ 1 ). (10) 
i=i 

An important step toward the solution of integrable models, in the thermodynamic limit, 
is the definition of new variables X(xj) = Xj, in terms of which S\l{xi, Xj) becomes a function 
only of the difference Aj — Xj. The corresponding integral equation derived from (9) will 
then have difference kernels. Following Baxter [25], we introduce a function 

w x 1, 1 + e~ 2r $(xi, x 2 ) s „ 12/ \ / 1 1 \ 

X(x 1 ,x 2 ) = - In 1 + e 2r$( Xi X2 } » ^i>^2) = 5' 12 (a;i,a; 2 ); (11) 

where r is the Baxter parameter, conveniently chosen for our models 

{et\ + J 2 cosh 2?7, for model A 
(12) 
cosh 2?7, for model B 

It follows (see [25]) the function X{x\,x 2 ) have the nice property 
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X(xi,x 3 ) = X(xi, x 2 ) + X(x 2 , x 3 ), (13) 

which imply 

\(xi,x 2 ) = \(xi)-\(x 2 ). (14) 

Using (11) and (14) we rewrite the Bethe-ansatz equations in the difference form difference 
form 



(jf)L=n ^-Ag>-r) ] = 1 ^ 
a =i smh(Aj - Xa +r) 



ft sinh(V-A^-r) " .mh(Ag»-Ag'-2r) 

J-l • w\ \(!) , a • umW \(i) , o n' i-,-,m. [10) 

j= i smh(A i — A« + r) sinh (A], ; - AV + 2r) 

where a and /3 are the arbitrary values. For example, we may choose a = and /3 = r for 

our convenience. The function $(x, 0) has the same form for both models, namely 

A(x ) = -ln — 7= / + r, $(rr,0 ) = — r- \ - 1 . 16 

The inversion of (16) gives us 



—J 1 cosh A sinh r ± y sinh A cosh r + cosh A sinh r tan 2 $ 

x = — • r~7T . x ; £ = 

smh(A + r) 



J 1 sinh A cosh r ± \ cosh 2 A sinh 2 r + sinh 2 A cosh 2 r tan 2 1? . . 

x= * — — r ; e = -l. (17) 

smh(A + r) v 7 

It is clear from (16) that the Bethe-ansatz equations have the . same form for . both models 
at the same values of the parameters r and . 

Let us consider the Bethe-ansatz equations in some limiting cases. At cos-# 1 we 
obtain Xj = sinh(Aj — r)/ sinh(Aj + r) and (15) gives us the Bethe-ansatz equations of the 
anisotropic supersymmetric t-J model, with anisotropy r [12]. In our derivation of (15) we 
assume the amplitudes in the eigenfunctions, corresponding to double site occupations, are 
related with those with single occupancy. Strictly at cos-# = 1, this assumption is not valid, 
unless there is no double occupancy as in the t-J model, and we should restrict m < L/2 in 
(15). 



In the limiting case of model B with £ = 1, where cos-# — > 0, rj — > 0, with [/ = 4?7 2 / cos-# 
fixed we obtain from (5) the Hubbard model with on-site interaction U — U. The relation 



(12) gives us r = cos(t?)/2 and by choosing Aj = z(7r/2 — 2 sin k yj cos($)/i7), A^ = 
i(it/2 — 2Aj^Jcos($)/U) we obtain the Bethe-ansatz equations of the Hubbard model [1] 

ik . L sinfc,- - A,- -iU/A . 1 , 10 . 

sin fcj- — hj + iU/A V ' 

"anhj-An + iU/A JL A p - A a + iU/2 

= — — = — ? = l,...,m. 19 

/i sin ^-^-$17/4 ^=1 A^ — A Q — iU /2 

The Hubbard limit can also be obtained in the limiting case of model A with e — 1 where 
cos-# — > 0,77 — > 0, but U — V 2 i — e 2ri cos($)/2 kept fixed. In this case we see from (5) that 
shifting Cj t 2 — * Cj-1,2, we recover the Hubbard model with on-site interactions U = V21. The 
Bethe-ansatz equations (18) are obtained from (15) by choosing Aj = i(7r/2 — 2e~ v sin fcj) 
and A« =i(7r/2-2e- T 'A a ). 

It is also interesting to observe that rational Bethe-ansatz equations can also be obtained 
for both models in the limit where r — > or r — > i7r/2. We should remark that even in this 
case we obtain new integrable quantum chains. For example at r — > we can rewrite (15) 
as 



e 



ikjL 



where 



II 7_7t! . i = l,-,n, (20) 

a=l A j lv a 2 
m \_A_i n A _ A _ *' 

TT ? A ? =-IT * 7 . « = l,...,m; (21) 



Aj = + 1) cot(%/2) + (J- 1 - 1) tan(%/2)], for e = +1, (22) 



and 



\j = [(J' 1 + 1) tan(%/2) + (J- 1 - 1) cot(%/2)]-\ for e = -1. (23) 

These solutions correspond to the model A at rj = 0,e = +1 and at cos d cosh 77 = ±l,e = 
— 1, and to model B at 77 = for both signs of e. 



To summarize we have presented two new two-parameter integrable models that general- 
ize the Hubbard and supersymmetric t-J models, and derived their Bethe-ansatz equations 
through the coordinate Bethe-ansatz method. Our results certainly motivate subsequent 
studies. One of them is the calculation of the phase diagram and critical exponents for 
arbitrary values of rj and Another interesting point raised by the present work, is the 
possible existence of a generalized R-matrix that reproduces that of the Hubbard model [27] 
at special points. It will be also worthwhile to generalize the model (1) for the case a > 2 
and to construct, in such way, the quite interesting Hamiltonian of multy-color Hubbard 
model [16,28]. 

We thank M. J. Martins and V. Rittenberg for useful discussions and also A. Lima-Santos 
and J. de Luca for interesting conversations. This work was supported in part by Conselho 
Nacional de Desenvolvimento Cientffico e Tecnologico - CNPq - Brazil. 
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